Thermodynamics of the Spin Luttinger Liquid in a Model Ladder Material 
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The phase diagram in temperature and magnetic field of the metal-organic, two-leg, spin-ladder 
compound (C5Hi2N)2CuBr4 is studied by measurements of the specific heat and the magnetocaloric 
effect. We demonstrate the presence of an extended spin Luttinger-liquid phase between two field- 
induced quantum critical points and over a broad range of temperature. Based on an ideal spin- 
ladder Hamiltonian, comprehensive numerical modelling of the ladder specific heat yields excellent 
quantitative agreement with the experimental data across the complete phase diagram. 

PACS numbers: 75.10.Jm; 75.40.Cx; 75.40.Mg 



Quantum spin systems display a remarkable diversity 
of fascinating physical behavior. This is especially true 
for systems such as spin ladders, which have a gapped 
or a gapless ground state, respectively, for an even or 
an odd number of ladder legs [1,. For two-leg ladders, 
and in general for any even leg number, quantum phase 
transitions (QPTs) between gapped and gapless phases 
can be driven by an external magnetic field. While these 
QPTs are generic in quantum magnets [2], the nature of 
the gapless phase depends crucially on the dimensional- 
ity of the spin system. In two and higher dimensions, 
a quantum critical point (QCP) separates the low-field, 
quantum disordered (QD) phase, with gapped triplet ex- 
citations, from a gapless phase with long-range antifer- 
romagnetic (AF) order, which can be well described as a 
Bose-Einstein Condensate (BEC) of magnons [HISJIl]. 

By contrast, for one-dimensional (ID) systems such as 
ladders, both long-ranged magnetic order and BEC are 
precluded due to phase fluctuations. In addition, spin ex- 
citations are best viewed as interacting fermions, whereas 
a bosonic representation pertains in higher dimensions. 
The physics of the gapless phase in ID is thus quite dif- 
ferent. It is a (spin) Luttinger liquid (LL) fS], and is a key 
component of the rich phase diagram presented in Fig. 1 
[31 El [II [3 [9] . In the LL, the spectrum is gapless with 
algebraically decaying spin correlations. Because there 
is no finite order parameter, the LL regime is reached 
from the high-temperature, classical regime through a 
crossover rather than a phase transition. Nevertheless, 
clear manifestations of LL behavior are expected not only 
in the correlation functions but also in thermodynamic 
quantities such as the magnetization and specific heat. 
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FIG. 1: Field-temperature phase diagram of the spin- 
ladder compound (Hpip)2CuBr4, showing quantum disor- 
dered (QD), quantum critical (QC), and spin Luttinger liquid 
(LL) phases. QCPs occur at Be (closing of spin triplet gap A) 
and Ba (spin system fully polarized). The contour plot shows 
the magnetic specific heat as Cm {T, B) /T. Local maxima 
from the reduction of the triplet gap by the Zeeman effect are 
indicated by crosses. Circles denote the LL crossover based 
on measurements of the magnetocaloric effect [Fig. 4], black 
lines are fits to extract the critical fields, and the dashed blue 
line indicates the onset of long-ranged order below 100 mK 
[211 122). Inset: lattice structure of (Hpip)2CuBr4 in projec- 
tion along the fo-axis, with Cu atoms blue and Br white. 



However, materials in which to explore such effects are 
rather rare. Investigations of the spin excitations and 
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thermodynamic properties of ladder compounds have to 
date been performed primarily on copper oxides. In these 
materials, the exchange interactions are typically some 
hundreds of meV, and thus the phases are not easily 
controlled by a magnetic field. Candidate ID materials 
in which thermodynamic properties have been studied 
around the QPT include the bond-alternating S — 1/2 
and S ^ I chains F5PNN [10] and NTEMP [11], the 
S = 1 Haldane system NDMAP [UlIS], and the 5 = 1/2 
system CuHpCl [14j , which was for some time considered 
to be a spin ladder |15j . While measurements in these 
materials show indications of LL behavior in parts of the 
field-temperature phase diagram, their magnetic proper- 
ties are influenced in large part by the presence of sig- 
nificant single-ion anisotropy and/or three-dimensional 
(3D) interactions between the chains, which tend to dom- 
inate the low-temperature specific heat at all fields [IB] . 
Hence we have pursued the search for materials suitable 
to study the intrinsic spin LL physics by seeking those 
showing, at minimum, a clear separation of energy (tem- 
perature) scales between ID and 3D interactions. 

Here we present the results of thermodynamic 
measurements on an exceptional two-leg ladder ma- 
terial, the compound piperidinium copper bromide 
(C5Hi2N)2CuBr4 [T71 ESI UHl US HD IMj, where all of 
the phases of interest can be accessed, as summarized 
in Fig. 1. In particular, we find an extended region in 
the phase diagram, between 0.1 K and 1.5 K, where a 
spin LL is observed. We show that the crossover into the 
LL is signaled by clear features in both the specific heat 
and the magnetization. In the gapless spin LL, the mag- 
netic specific heat is linear at low temperatures. Its field- 
and temperature-dependence are in excellent agreement 
with numerical calculations involving no free parameters. 
This demonstrates that the material is very accurately 
described by a minimal spin-ladder Hamiltonian. 

High-quality single crystals of (C5Hi2N)2CuBr4, ab- 
breviated (Hpip)2CuBr4 in the following, were grown 
from solution. In this material, the 5=1/2 magnetic 
moments of the Cu^"*" ions are arranged in a ladder-like 
structure along the a-axis [Fig. 1, inset]. The rungs (Jr) 
of this ladder are formed by two equivalent Cu-Br-Br- 
Cu superexchange paths with a center of inversion sym- 
metry [in], while the legs (J;) involve one similar but 
longer interaction path. The ladder units (Cu2Br8)'^~ are 
well separated by the organic (C5Hi2N)+ cations, which 
contribute only very little to the electronic properties of 
the host structure, and hence any magnetic exchange be- 
tween ladders (J') is expected to be small. Direct mea- 
surements of these interactions, Jr=12.9(2) K, J;=3.3(3) 
K, and J' < 100 mK, based on inelastic neutron scatter- 
ing experiments [23 , are in very close agreement with the 
values extracted from magnetostriction [20j and nuclear 
magnetic resonance (NMR) [21^ measurements. 

The specific heat and magnetocaloric effect (MCE) 
were measured on a purpose-built calorimeter at the 
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FIG. 2: (a) Measured total specific heat CtotiT) at fixed mag- 
netic field. Solid (dashed) lines are based on fits to the data 
before (after) background subtraction. The blue line indicates 
the uniform, non-magnetic background, (b) Magnetic specific 
heat CmiT,B) for B < B^, and (c) Be < B < B^. Lines in 
(a-c) show , and are based on ED and DMRG calcu- 
lations, as explained in the text, (d) Cm{T, B) /T measured 
at fixed temperature. The region with linear temperature- 
dependence of the specific heat is indicated by LL, while T 
marks the peak due to the softening triplet. 



Helmholtz Centre Berlin (Laboratory for Magnetic Mea- 
surements at BENSC), using single crystals of masses 
4.78 mg and 9.69 mg in the respective temperature and 
field ranges 0.3 K to 15 K and T to 14.5 T. The field 
was apphed parallel to the crystallographic a-axis, a ge- 
ometry in which we obtained the values Be = 6.99(5) 
T and Bs = 14.4(1) T for the two QCPs [Fig. 1]. The 
specific heat was extracted from a quasi-adiabatic relax- 
ation technique and, using the same set-up, the MCE 
was recorded with a sweep rate of 0.05 T per minute. 

In Fig. 1 the magnetic component of the specific heat 
Cm/T of (Hpip)2CuBr4 is presented across the entire 
phase diagram. It shows clearly three distinct regimes: 
QD, quantum critical (QC), and spin LL. The contour 
plot was obtained from 27 scans in field and temperature, 
after subtraction of a field-independent lattice contribu- 
tion C/(T) and of a small nuclear term, which both are 
determined from a simultaneous fit to all available data. 

In Fig. 2 (a)-(d) we show individual measurements 
of the total and the magnetic specific heat, respectively 
Ctot{T) and C^{T,B). In the QD regime, B < B^, C„ 
shows a single peak at approximately 5 K [Fig. |2jb)]. 
This peak is attributable to the triplet excitations of the 
ladder, and is exponentially activated at lower tempera- 
tures due to the presence of the spin gap A [24] . With 
increasing field, the gap is reduced by the Zeeman effect 
(A — > A — gjiBB). Field scans such as those in Fig. 2(d) 
show most clearly the reduction of the gap and are used 
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to extract the critical field Be in Fig. 1, yielding very 
good agreement with determinations by complementary 
experimental techniques p i [19 1 [20 l [2T| [22] . 

The specific heat changes dramatically for fields B > 
Be, which we explain by the formation of the LL phase 
[Fig. 2(c)]. While at high temperature Cm is domi- 
nated by the (gapped) Sz — triplet states, at low 
temperature an additional peak develops. Below this 
peak, the temperature-dependence remains linear up to 
B — Bg K, 14.5 T, with a field-dependent slope. The 
linearity of Cm is demonstrated in Fig. |2jd) . For fields 
near the maximum of the LL dome, the ratio Cm IT mea- 
sured at different temperatures collapses onto the same 
curve. This temperature-dependence is consistent with 
the presence of gap less spinon excitations with a finite 
velocity m, its slope being inversely proportional to u 
[5]. The first peak thus occurs when the temperature is 
large enough to probe deviations from this linear regime. 
It can be taken as an estimate of the crossover to en- 
ter the LL, and is visible in Fig. [l] for Be < B < Bs- 
The field-dependence of C,„ is almost symmetric about 
{Be + Bs)/2 « 10.7 T. In the strong-coupling limit, 
Jr/Ji ^ 1, a perfect symmetry would be expected due to 
the exact particle-hole symmetry of the XXZ chain in a 
field [Mj [25] . Here we observe clearly deviations charac- 
teristic of the underlying ladder structure. Similar effects 
are also visible in spin correlation functions and in the 
low-temperature phase diagram, which can be measured 
by NMR gl] and neutron scattering ^ . 

At B > Bs, the specific heat becomes exponentially ac- 
tivated again due to the opening of a field-dependent spin 
gap in the fully saturated phase. However, this regime 
is close to the limit of our experimental window, and so 
the high-field phase is not investigated further here. 

The experimental data have been compared with 
several theoretical calculations, and the agreement is 
remarkable (Fig. 2). Numerical results were ob- 
tained by exact diagonalization (ED) and by adaptive, 
time-dependent density-matrix renormalization-group 
(DMRG) calculations [25], both performed for a single 
ladder with Jr=13 K, Jr/Ji=^, and g=2.Q& (i.e. no free 
parameters). We stress that in both techniques it is im- 
portant to retain a sufficient number of ladder states for 
a quantitative description of thermodynamic data. The 
DMRG calculations (2x40 spins) may be regarded as the 
definitive behavior of this model. In the ED calculations, 
the specific heat of even- (odd-)length ladders converges 
rapidly from above (below) to the infinite-size limit; thus 
finite-size effects are essentially removed here by taking 
an average between systems of 2 x 10 and 2x11 spins. The 
ED and DMRG results are indistinguishable both in the 
QD phase [Fig. 2(b)] and in the LL regime [Fig. 2(c)]. 
Slight deviations from the experimental data are found 
only close to the upper critical field Bg and at IIT. 

Some physical insight into the numerical results is 
afforded by two approximate treatments. A statisti- 
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FIG. 3; Comparison of calculated ladder specific heat, (a) ED 
vs. TTW for magnetic-field values B < Be [Fig. 2(b)]. (b) 
ED vs. BA for some magnetic-field values B > Be [Fig. 2(c)]. 



cal Ansatz (TTW) [24] developed for spin ladders, and 
shown previously to describe very accurately the thermal 
renormalization of triplet excitations in the 3D dimer sys- 
tem TlCuCla [57j, uses the correct number of states but 
applies their hard-core constraint only globally. In the 
spin ladder, this approximation underestimates the local 
energy of the excited states, leading to a systematic shift 
of weight to lower energies as B ^ Be [Fig. [3ja)]. A 
mapping of the lowest two modes of the ladder Hamil- 
tonian onto an effective 5=1/2 XXZ chain O [5S], 
whence thermodynamic quantities are computed exactly 
from the Bcthc Ansatz (BA), is very accurate for the 
low-energy physics ai B > Be, but cannot reproduce the 
heat capacity at higher temperatures because of the miss- 
ing triplet states (5*2 = 0,-1) [Fig. [3|b)]. We conclude 
that the thermodynamic properties of (Hpip)2CuBr4 are 
described very accurately by a model of a single two- 
leg ladder, and that comprehensive measurements of the 
specific heat identify an extended LL regime. 

We turn now to a different observable, the uniform 
magnetization (M), which is notoriously difficult to mea- 
sure at temperatures below 1.5 K. Very precise measure- 
ments can be obtained by NMR [21], but here we use 
an alternative method to probe the crossover to the LL. 
We determine the derivative of the magnetization with 
respect to temperature using the relation {6Q/SB)/T — 
— {dM/dT)\B, where SQ is the amount of heat created or 
absorbed by the sample for a field change SB due to the 
MCE. Figure |4] shows both representative {SQ/SB)/T- 
traces (corrected for a small base-line drift at higher 
temperatures) and a contour plot of all available data, 
presenting directly dM/dT. In the free-fermion model, 
which is an excellent qualitative description of spins near 
the QCP in ID [^, and in more refined approaches 
[S] [21122], the magnetization presents a minimum or max- 
imum as a function of temperature {dM/dT = 0). These 
extrema occur when the temperature matches the chem- 
ical potential, and thus provide another determination 
of the crossover temperature for the LL phase. The ex- 
tracted phase boundary and the positions of the peaks 
in the specific heat agree well within expectations, as 
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FIG. 4: Magnetocaloric effect in (Hpip)2CuBr4. (a) Heat- 
flow SQ to and from tfie sample as a function of magnetic 
field divided by temperature, (SQ/SB) /T. (b) Contour plot 
of dM/dT as a function of field and temperature. White 
circles denote the phase boundary derived from dM/dT — 
(see also Fig. 1), while black circles are maxima in the specific 
heat, dCm/dT = 0, obtained at fixed field. 

demonstrated by the solid symbols in Fig. |4|b). 

The structures present here in the magnetization and 
specific heat differ marfcedly from those occurring when 
there is a BEC. In that case, because a real phase tran- 
sition occurs, the specific heat diverges with a A-type 
anomaly. Such a shape has been observed at BEC transi- 
tions in higher-dimensional materials such as BaCuSi206 
[28j . At the same temperature the magnetization devel- 
ops a minimum, but with a cusp-like structure , which 
has been observed in TlCuCls g]. For ID spin ladders, 
the magnetization minimum and the specific-heat peak 
have a different origin: they correspond to the crossover 
to the LL regime. Thus there is no divergence in the spe- 
cific heat and the magnetization minimum is analytic, 
reflecting the absence of a phase transition; the tempera- 
tures of the two features, although similar, are not iden- 
tical. For (Hpip)2CuBr4, a real phase transition of BEC 
type does occur at a much lower temperature, Tn k 100 
mK (Fig. 1, [211122), due to a 3D coupHng of the ladders. 

In summary, we have measured the specific heat 
and magnetocaloric effect in the metal-organic, two- 
leg spin ladder (C5Hi2N)2CuBr4. The excellent low- 
dimensionality and optimal energy scale of the exchange 
interactions make this material unique, and allow a de- 
tailed investigation of the phase diagram in temperature 
and in fields up to magnetic saturation for the quantum 
spin ladder. We find an extended region of spin Luttinger 
liquid behavior over at least one order of magnitude in 
temperature, lying clearly above any three-dimensional 
physics triggered by residual inter ladder interactions. 
The high-precision experimental data have been ana- 
lyzed using the most advanced exact diagonalization and 
density matrix renormalization group techniques to cal- 
culate thermodynamic quantities for all of the phases 
(i.e. across two quantum critical points). From the di- 
rect and parameter-free fit of the experimental and nu- 
merical results, we conclude that (C5Hi2N)2CuBr4 is re- 



markably well described by a minimal spin-ladder Hamil- 
tonian, with other possible effects (frustrated interac- 
tions, Dzyaloshinskii-Moriya terms, lattice coupling) be- 
ing very small. Hence the material offers unprecedented 
opportunities to investigate the intrinsic physics of low- 
dimensional quantum systems. 
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